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Calculations for an inhomogeneous rf-wave superconductor 
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We calculate the distributions P[Ao(r<)] of local d- wave pairing amplitude Ao(rj) at a position 
inside a disordered high temperature superconductor (HTS) family. To reproduce the observed inho- 
mogeneity a random potential V lmp , within a Bogoliubov-de Gennes (BdG) formalism, is considered. 
We perform calculations with fixed values of the disorder strength V lmp obtaining the distribution 
of Ao(rj), and local density of charge carriers p(ri), for different compounds of the LSCO family. 
The calculation of the relative root-mean-square deviation shows that the underdoped compounds 
are more inhomogeneous than the overdoped ones. Also, the spatial variation of Ao(ri) decreases 
as the average density of charge carriers (p) increases, demonstrating that the compounds are more 
homogeneous for high values of (p) . The results indicate that the d-wave superconducting gaps seem 
to be more sensitive to a change in the disorder in comparison with the s-wave superconducting 
gaps. 

PACS numbers: 74.72.-h, 74.20.-z, 74.80.-g 



I. INTRODUCTION 

It is clear by now that HTS have many non- 
conventional physical properties, and they are the reason 
why the fundamental interaction is yet to be discovered. 
In some cases it is possible that these unusual properties 
result from the variation in the average density of holes 
(or charge carriers) (p) in the CuC>2 planesi^. 

Many experimental features show that the density 
of charge carriers (p) plays an important role in the 
physics of the HTS. It is observed that the zero tem- 
perature superconducting gap Ao increases when (p) di- 
minishes^&i, which is an unexpected behavior since the 
critical temperature T c vanishes for low values of 
It is also experimentally observed that for some com- 
pounds the electrical charges are highly inhomogeneous 
in the CuC>2 planes 6 -' 9 -' 1 ^. These intrinsic charge inhomo- 
geneitiesii are not related to crystal defects. Although 
it may be not so strong in some cuprate superconduc- 
tor o 12 : 13 ' 14 : 15 . The origin of the inhomogeneity observed 
in these materials may be related to the fact that the 
holes injected into the CuC>2 planes by chemical dop- 
ing may leave behind charged dopant ions 1 ^. Such in- 
homogeneity may produce regions with spatially vary- 
ing Ao(rj) and a variation in the local density of states 



(LDOS) 1 ?. Many similar results from scanning tunneling 
microscopy were reported lately on clean surface a 18 ' 19 . 
Inside the bulk, Nuclear quadrupolar resonance (NQR) 
experiments 20 have also measured that the charge inho- 
mogeneity increases as the temperature decreases, ex- 
actly as one expects in a phase separation transitior 1 '" 



Since an inhomogeneous medium does not have trans- 
lational symmetry, we applied the Bogoliubov-deGennes 
(BdG) formalism 3 - ! 16 : 21 : 22 : 23 : 24 , conceived originally to 
deal with finite systems. To reproduce the intrinsic 
charge inhomogeneity we consider a nonmagnetic local 
disorder potential, which has only the effect of changing 
the local chemical potential on each site of the lattice^ 1 .. 
Our purpose here is to obtain a distribution of local su- 
perconducting gaps Ao(rj), and a distribution of local 
density of charge carriers p(vi) for an entire HTS family. 
We show that with the same magnitude of the disor- 
der strength, the underdoped compounds become more 
inhomogeneous than the overdoped ones, which is in ac- 
cordance with neutron diffraction experiments^. For the 
(i-wave gap we observe that for very low disorder strength 
the distribution of local superconducting gaps gains sig- 
nificant weight near Ao(ri)=0 and the distribution of lo- 
cal charge carriers starts to bifurcate. This behavior dif- 
fers from the s-wave gap symmetry distributions 21 , which 
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needs higher values of disorder strength to have similar 
behavior. 



II. THE HAMILTONIAN 

To describe the charge carriers dynamics in the CuC>2 
planes of the HTS we consider an extended Hubbard 
Hamiltonian in a square lattice 

H = ~ tij4 a c ja + ^2(V- mp - \i)n ia 

{ij}a icr 

+U^2niinn + — Y n i<T n i<r', (!) 

i (ij)a-a-' 

where cl a (ci a ) is the usual fermionic creation (annihi- 
lation) operator at site r^, with lattice parameter a—1 
and spin o"{jJ.}. n ia — cj^c^ is the density operator, 
and Uj is the hopping between sites i and j. U is the 
magnitude of the on-site repulsion, and V is the magni- 
tude of the nearest-neighbour attractive interaction, fx 
is the chemical potential, and V- mp is the magnitude of 
the local disorder potential at site i, defined here as a 
random independent variable, uniformly distributed over 
[-V imp , V imp ], similar to what Ghosal et alSL made to a 
S-wave superconductor. The value of V- mp assigned to 
the site controls the strength of the disorder. We do 
not know the origin of V* mp , but it is possible to spec- 
ulate its origin: It can be due to phase segregation, as 
in the case of I^CuC^+a, which has been observed to 
form oxygen rich and poor phase o 25 ' 26 , or ion diffusion as 
observed in microcrystals of Bi22132i. These findings, to- 
gether with the NQR experiment of Singer et al2£, led us 
argue that is very likely that the upper pseudogap line is 
the onset of phase separation or bimodal decomposition 3 
in Cuprates superconductors. However, independently of 
its origin, the effect of V i enters in our calculations by 
the change of the local chemical potential at each site i 
and, as a consequence, the local density of charge carriers 
p(r.;), gives rises to the charge inhomogeneity. Whether 
such charge non-uniformity is in form of stripes 2 ^ is still 
an open question. There are many ways of introducing 
the effects of disorder, for instance, Ovchinnikov et al22 
considered magnetic impurities and Nunner et ali^ as- 
sumed a disorder in the pairing strength coupled with 

t Amp 



III. THE LOCAL EQUATIONS 

To apply the BdG theory to the Hamiltonian (fl]), one 
may define the pairing amplitudes, or pair potentials 23 

A s (Ti) = V{cnci + s<t), &u(ri) = C/(c a c iT ), (2) 

where 5=±x,±y are unit vectors for a square lattice. In 
the mean-field theory Eq.fT]) can be solved using these 



pairing amplitudes. The resulting effective Hamiltonian 
H e ff is given by 

iScr icr 

+ X)[ A «( ri ) Ci ± Ci +' 5 t + A *( r i) c ]+5T c li] 
iS 

+ ^[Ac;(r,)4 T 4 4 + AUri)^}, (3) 

i 

where fMi=Li-U/2(rii) is the local chemical potential, 
which incorporates the site dependent Hartree shift 
■j(rii)2i. Both, the Hartree shift and y" np take care of 
the charge inhomogeneity of the system. The electronic 
density is given by (ni)=E cr (rii cr }, and the hole density is 
p(ri)=l-(m). The hole type density of charge carriers of 
a specific compound is given by 

1 Ns 

i—l 

where N s is the number of sites of the N x N square 
lattice. Therefore, Eq.(|3]) fixes the chemical potential /i. 
The H e ff is diagonalized by the BdG transformations 21 

C A = E„[7nT"n(ri) -7«j.<( r *)]i 

C H = ErJ7nl««( r i) +7l T <(r*)], (5) 

where -f n(T and 7^ are quasiparticle operators. u n (rj) 
and v n (ri) are normalized amplitudes for each r^, and 
are obtained from the BdG equation a 21 i 23 

(K A \ /u„(r,)\ fu n (ri)\ 

\A* ~K*J \VniTi)) n \VniVi)) 

with 

Ku n (vi) = -^2t l . l+s u n (r l + S) + iyl mp - L~Li)u n {Ti) 

s 

Au n (vi) = ^2 A s {r l )u n (r l + S) + A c/ (r l )u n (r l ), (7) 

<5 

and similar equations for v n (r^) . These equations give the 
quasiparticle eigenenergies E n (> 0). The temperature 
dependent BdG equations are solved self-consistently to- 
gether with the pairing amplitude o 23 i 24 

A u (r i ) = -l/EnttnCriKfrOtaiihgj^, (8) 

+<(r J )u„(r. i + <S)]tanh 5 f^ T . (9) 

Eq.© has four different possibilities of directions for a 
square lattice: A±x(ri) and A±y(r.i). The combination 
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o Ino et. al (ARPES) 
■■■■ V imp =0.1t (min. gap) 
- V imp =0.0 

- V imp =0.1t(max. gap) 



20 - 




<P> 



FIG. 1: The experimental leading edge gap of Ino et al.— (cir- 
cles), interpreted as a superconducting gap, and the Ao curve 
as a function of the average density (p) for the L^-^Sr^CuCU 
family with d-wave symmetry. The solid line is the result 
without disorder (V mp —0). The long-dashed and doted lines 
correspond to maximum and minimum pairing amplitudes, 
respectively, for the disorder strength V tmp —0.1t. 
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FIG. 2: The distribution P[Aofa)] of local d-wave super- 
conducting gaps Ao(rj) for the near optimum compound 
(p}=0.15, for different values of the disorder strength y tmp . 
At low disorder, P[A(rj)] is peaked around the average gap 
value Ao=9meV. As the disorder increases, P[A(rj)] becomes 
broad. At very large disorder P[A(rj)] gains significant weight 
near Ao(i"i)=0. 



of these terms may give rise to a d-wave gap or s-wave 
gap 23 . The d-wave case is given by 

A fa) = ±[A*fa) + A_ x fa) - A y fa) - A_yfa)]. 

(10) 

Also, the local hole density of charge carriers is given by 
p(r,) = l-2^[| Mn (rO| 2 /n + |^(r l )| 2 (l-/n)] ! 

n 

(ii) 

where f n is the Fermi function. With Eq. (fTT|) one can cal- 
culate pfa) which, together with Eq.flU), gives the density 
of charge carriers (p) of a compound. 

Therefore, the BdG equations are solved self con- 
sistently, together with Eq.© for Ay fa), Eg. ([10]) for 
Aofa), and with Eq. irpTj) for pfa), with periodic bound- 
ary conditions on a lattice with N s sites. We have per- 
formed calculations with lattices from 14 x 14 to 24 x 24, 
but here we concentrate on the parameters that repro- 
duce the experimental results in a 22 x 22. 

IV. THE LOCAL DISTRIBUTIONS 

In order to reproduce the experimental phase dia- 
gram' 51 Aq us (p) of LSCO, we considered values of U 



and V close to values used in literatu re 23 ' 2 . The con- 
stant coupling U has low influence in the calculations 
of the superconducting gap for the d-wave symmetry^. 
The hopping integrals were based on ARPES experi- 
ments results 3 ^ and we considered hopping integrals from 
the first to the fourth neighbour, differently from previ- 
ous works with the BdG formalis m 3 ' 16 ' 21 ! 22 ' 23 ' 24 . Here, 
the average superconducting gap is written as Ao, dif- 
ferent from the local superconducting gap, Aofa). The 
Aofa) gap can be also obtained by a different procedure 
through the study of the local density of states (LDOS) 34 , 
which is given by N t {E) = £ n [K(x;)| 2 ' f' n {E - E n ) + 
\v n (xi)\ 2 f n (E + E n )], where the prime is the derivative 
with respect to the argument and the E n are calculated 
by the BdG equations (Eqj6]). These local calculations 
are important to interpret the new STM results which 
show differences in the LDOS at mesoscopic scal o 17 ' 18 ' 19 . 

In FigQ]we plot the results for the average d- wave gap 
Ao, for a 22 x 22 lattice with the values: <=ii=0.35eV, 
t 2 =0.50£, i 3 =0.27t, f 4 =0.15i, £/=0.35i and V=-0.45t. 
These values are very close to previous calculations and 
known to reproduce the curve of T * 14 ' 15 . From FigQ] we 
observe that Ao diminishes as the density of charge car- 
riers (p) increases, which is a common behavior of the 
superconducting gap for most HTS^i 7 -' 3 -^. As mentioned 
T* vs (p) 7 - ' 14 ' 15 decreases from high values for low densi- 
ties (underdoped region), into low values for high densi- 
ties (overdoped region). Since experimental results 5 -! 7 -^ 
and theoretical calculations 14 ' 15 indicate that the super- 
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FIG. 3: The distribution P[p(ri)] of local density of charge 
for the near optimum compound, with average density value 
(p}=0.15, for different values of the disorder strength y tmp . 
For low values of V lmp , P[p(ri)] is peaked around the average 
value. For V imp =0.0A5t and 0.070t, P[p(r=.)] seems to have 
a bifurcation. For V lmp greater than O.li, P[p(r 4 )] tends to 
spread. 



As we have already discussed in Refs . 14 i 15 , the maximum 
local gap is measured in the tunneling experiments and 
it seems to increase with the level of inhomogeneity 4 . 
From Figf3]we see that, indeed, the disorder spreads the 
distribution P[A(r=.)], what increases the maximum lo- 
cal superconducting gap. Since the experimental leading 
edge shift, interpreted as the magnitude of the supercon- 
ducting gap 34 , increases dramatically in the underdoped 
region, we can speculate that this is possibly the result of 
the inhomogeneity in the Cu02 planes of the underdoped 
compounds. One can see the effect of V unp =0.lt on the 
maximum local superconducting gap curve in FigfT] 

It is important to notice that, since the values of y™ riv 
are chosen randomly at each site between the values in 
the interval [—V mp , V lmp ], the results are averaged over 
10 realizations of disorder. With this procedure we have 
reproduced the results of the s-wave gap in a two dimen- 
sional 12x12 lattice of Ref^. 

In Fig|3]we plot P[p(ri)) for the (p)=0.15. We observe 
that, at low disorder, P[p(r,)] is peaked around the av- 
erage density (p)=0.15, and for V imp =0M5t and 0.070i, 
P[p(ri)] seems to have a bifurcation. As the disorder 
increases, P[p(vi)\ also spreads considerably. 



V. THE INHOMOGENEOUS MEDIUM 



0.2 



conducting gap Ao of the HTS scales with T* and not 
with T c , it is reasonable to expect a similar behavior of 
T* vs (p) and Ao vs (p). The oscillations of the the- 
oretical curves in FigfT] are associated with the number 
of hopping integrals: To embrace the entire range where 
the superconductivity of LSCO exists, it was necessary 
to consider hopping integrals from the first to the fourth 
neighbour, which generates the oscillations in the Ao vs 
(p) phase diagram. 

To study the disorder effects, with the same coupling 
parameters and hopping integrals used in the homoge- 
neous case, we calculate the distributions of local d-wave 
superconducting gaps P[Ao(i"i)] and local hole densities 
P[p(ri)]. In FigEJwe plot P[A (r. t )] for the (p)=0.15 
compound. We observe that, for low disorder, P[A (rj)] 
has a sharp peak near Ao(i"i)«9meV, which is the homo- 
geneous result, in accordance with Fig[JJ As the disorder 
increases up to l.Oi, P[Ao(rj)] becomes broad, and the 
system becomes more inhomogeneous. At very large dis- 
order P[A(r,-)] gains significant weight near Ao(ri)=0. 
This behavior is similar to the s-wave gap, although it 
is necessary a disorder of 2t in a 12x12 lattice, and 3t in 
a 24x24 lattice, to obtain the same result for the s-wave 
cas o 21 ' 36 . Therefore, the d-wave gap is more sensible to 
changes in the disorder than the s-wave gap. As one ob- 
serves from FigfSJ for V™ mp =0.045i the disorder generates 
a distribution of Ao(r=,) around the average gap value Aq. 
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FIG. 4: The distribution P[Ao(n)] of local superconducting 
gaps for some selected compounds: the underdoped (/o)=0.10; 
the near optimum (p)=0.15, and three overdoped compounds, 
(p}=0.20, 0.25 and 0.30. The magnitude of the disorder is 
V mF =0.Wt. From the figure we observe that, as the aver- 
age density (p) increases, the distributions P[Ao(i"j)] become 
sharp and thin. 



In this section we study the inhomogeneous case of 
V imp =0.10t. 
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FIG. 5: The distribution P[p(n)] of local charge carriers 
for: the underdoped (p)=0.10, (p)=0.15, 0.25 and 0.30. 
V imp =0.10t. Note that as (p) increases, P[p(n)] becomes 
more peaked. 
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FIG. 6: The relative root mean-square deviation rms.rel for 
LSCO with V lmp =0.1t. It indicates that as (p) increase, 
rms.rel decreases showing that the local gap distributions be- 
comes narrower. 



In FigH] we plot the distribution P[Ao(r.j)] for some 
selected compounds for V lmp =Q.10t. From the figure we 
see that as the average density (p) increases, P[A (rj)] 
becomes sharp and thin. Therefore, it is clear that the 
average variation in Ag(rj) decreases with increasing av- 
erage doping. In FigOwe plot the distribution P[p(rj)] 
for V' Lmp —0.1Qt. Again we observe that increasing (p) 



the distributions become sharper. 

a) 





FIG. 7: Spatial variation of Ao(ri) in the 22x22 lattice, for 
the mean densities (p)=0.15 in a), 0.20 in b), and 0.30 in 
c). It was considered the same realization of disorder for all 
compounds. In the figures the light regions are high Ao(r,) 
values. Note that for low value of (p) there are many different 
color regions, indicating a more inhomogeneous medium. 



To study the range of inhomogeneity, we calcu- 
late the relative root-mean-square deviation (rms.rel), 
[((AA ) 2 )]V A o, where, by definition, [((AA ) 2 )]2 gives, 
approximately, a linear measure of the width of the 
range over which Ao is distributed. In order to com- 
pare the rms.rel of each LSCO compound, the distribu- 
tions P[Ao(r.;)] must be shifted, so they all have the same 
mean gap value A , maintaining their original shape. We 
can see in Fig|5]that, as the average density (p) increases, 
rms.rel diminishes, i.e., the range over which the super- 
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conducting gap Ao is distributed also diminishes, which 
indicates again that as (p) goes into the overdoped re- 
gion, the compounds become more homogeneous. 

In Fig[7|we show the spatial variation map of Ao(rj), 
in the 22x22 lattice, for three different compounds 
for a particular realization of the random potential 
(V lm P=0.10t): in FigJT^ we have the near-optimum com- 
pound ((/j)=0.15), and in Fig[7)3 and c, we have two over- 
doped compounds, (p}=0.20 and 0.30, respectively. In 
the figures the light regions corresponds to high local su- 
perconducting gaps, and the dark regions corresponds to 
low local superconducting gaps. In Fig(7£i we observe the 
coexistence of light regions and dark regions, indicating 
a high inhomogeneous regime. As the density increases, 
most of the regions become darker, indicating a more ho- 
mogeneous medium. Therefore, we conclude again from 
the figures that, as the density (p) increases, the systems 
tend to become more homogeneous, in agreement with 
FigJU What means that a d-wave superconductor be- 
comes more homogeneous as the concentration of holes 
or doping level increases. 

VI. CONCLUSION 

We have calculated distributions of local (i-wave su- 
perconducting gaps Ao(rj) and local density of charge 
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